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Abstract 



I We present a numerical study of the ghost propagators in 

Landau gauge for SU(2) and SU(3) gauge theories at (3=2.7 
. and /3=6.0, respectively. Analyzing different lattice sizes up 

to 32^, we find small finite size effects. Down to the small- 
est available momenta, we observe no evidence for dipole 
^ . behaviour of the ghost propagators. 

> 

X 



c5 ' 1 Introduction 



The gluon propagator, P, albeit a non-gauge-invariant quantity, is considered 
of prime interest in the quest to gain intuitive insight into the physics of con- 
finenent in non-Abelian gauge theories. From a linearly rising potential, one 
would argue that P is dominated by a singular ity|j|]. This infrared be- 
haviour can be drastically altered, however, in a scenario of dynamical gluon 
mass generation as pointed out in Ref.[^. 

The issue is closely related to the singularity structure of the Green function 
of the Faddeev- Popov ghost. In his seminal paper on gauge fixing ambiguities 
Gribov Q has dealt with the implications of the proper choice of integration 
range (in the functional gauge field integrals) onto the singularity structure of 
the ghost propagator. This is elaborated in detail in a recent comprehensive 

^We ackowledge A. Hulsebos who contributed substantially during the early stages of this 
work. 
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paper by Zwanziger |^ where the theoretical arguments are presented in quite 
some detail that lead up to the prediction of a dipole-type infrared singularity 
in the Faddeev-Popov propagator. 

It is well accepted by now that the lattice laboratory provides a valuable tool to 
the heuristic study of this highly intricate non-perturbative scenario. Indeed, 
lattice investigations of the gluon propagator in Landau gauge did reveal a mo- 
mentum dependent effective mass for the gluon|^, It should be remembered, 
though, that - in view of the very Gribov ambiguities - such lattice results rely 
deeply on the quality and efficiency of gauge fixing procedures on the lattice 




In this letter, we wish to complement the previous gluon propagator investiga- 
tions by presenting a first numerical study of the ghost propagator via direct 
lattice simulation. 



2 Faddeev-Popov Operator 

The ghost fields can be defined on the lattice^ in a similar manner as in the 
continuum. The ghost propagator G is given in terms of the inverse of the 
Faddeev-Popov operator M: 

M = -VD{U), (1) 

where D[U) stands for the covariant derivative on the lattice. It has a diagonal 
form in the algebra space, being an average over field configurations: 

G{x-y)5'^^^{M-^''\x,y)). (2) 

In a Landau gauge fixed configuration the action of the operator M on an 
arbitrary element w in the algebra space A of the gauge group SU(N) is given 

by 



ir^^ [A%x)u:\x + fi)- 4(x - i,)u,^ix - A) ] }. (3) 



2 

Here 5" is a linear functional 



S;%x) = --TV [(rV" + rV) (f/^(x) + U;{x)) J . (4) 

yl^ is the gluon field and the Pauli matrices r are used to span an antihermitian 
basis of the linear space A. 



*A very useful collection of the relevant lattice formula related to the ghost fields can be 
found in Ref. fll. 
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3 The Simulation 



Run Parameters. We ran simulations on a series of lattices, whose param- 
eters are summarized in Table |l[ 





size 


a 


# Configurations 


SU(2) 


16* 
32^ 


2.7 
2.7 


50 
15 




8^ 


6.0 


64 


SU(3) 


16^ X 32 
24^ 


6.0 
6.0 


30 
10 



Table 1: Lattices used in this simulation. 

The SU(2) gauge configurations were generated using the Kennedy-Pendleton 
heatbath algorithm. 2000 updates were discarded for thermalization, while 
gauge fixing and ghost propagator calculations took place every 100 updates. 
For SU(3), we applied a hybrid algorithm of Cabibbo-Marinari heatbath and 
Creutz overrelaxation steps with mixing probability 1:5. The configurations 
have been analyzed every 250 sweeps after 2000 thermalization steps. 

Gauge fixing. Gauge fixing was done by minimizing the functional 

^t/b] =EE ; Uf,{x)=g{x)U^{x)g\x + fL). (5) 

X IJ, ^ ' 

This lattice condition is slightly stronger than the conventional Landau condi- 
tion in continuum theory: 

d^^A^ = 0. (6) 

The minimization procedure may be carried out by use of one of the standard 
relaxation algorithms which will drive the system to one of its local minima of 
F;7[g], thus delivering one of the possible Gribov copies. 

We aim to achieve configurations lying in the fundamental modular region A, 
which is given by the set of absolute minima of Fu\g\ on all gauge orbits. 
The way to do that is to modify the gauging algorithm and to mobilize the 
system such as to allow - in a gentle manner - for escape from the attraction 
of its current closest minimum. To that end we applied two approaches: (a) 
the simulated annealing (SA) for SU(2) and (b) the stochastic overrelaxation 
(SOR) algorithms for SU(3): 

(a) The idea of the annealing algorithm is most easily exposed in the language 
of spin models: Minimizing -F[/[(7] may be viewed as retrieving the ground state 
of a spin system, with action being given by 

Sv{g) = -(5sFu[9\. (7) 
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In this picture, the ground state may be reached by performing Monte Carlo 
sweeps on the 'spin' degrees of freedom g. 

Let us give our prescripton of the anneahng algorithm in detail: (1) start at a 
spin model coupling (5s = 0.20, and perform 10 Creutz or Kennedy-Pendleton 
heatbath updates on the spin model, at given (3s\ (2) increase fig in steps of 
A/?s = 0.20 and keep updating until fig = 25.0 is reached; (3) at this point turn 
to ordinary relaxation. 

This gauging procedure was found to be very efficient on both 16^ and 32^ 
lattices. 

(b) In the stochastic overrelaxation, a "wrong" gauge transformation is applied 
once in a while, i.e. with probability w, during the iteration process. If w is 
sufficiently high, the system may be driven away from a local minimum. We 
applied this method in conjunction with the Los Alamos algorithm |8| after 
mapping the lattice onto a checker board basis. We have set the probability 
u;=0.7. Our criterion for achieving the Landau condition was 

la^^^^l^ < 10"^ (8) 

On our largest SU(3) lattice, 24^, about 1200-1400 iterations were needed. 



Inversion of M. In Landau gauge, M = — VZ? = —DV is a singular matrix. 
It annihilates all zero modes on the left and on the right hand sides. Therefore, 
in order to compute its inverse, we have to separate the zero modes and perform 
the inversion on its regular part only. M is a real symmetric matrix acting on 
the algebra space A of the gauge group, which may be decomposed into two 
subspaces: 

A = AoeAi (9) 

such that 

MAo = and M Ai ^ (=0 only for the null vector). (10) 

We start with the simple observation that the vector Mv lies in for arbitrary 
V and compute G by solving the linear equation 

MMv = Ms, (11) 

using the conjugate gradient (CG) and minimal residuum (MR) as standard 
iterative algorithms. Both methods start off from some initial guess vector vq 
and achieve an improved approximation to the solution by adding a vector dvi 
to the current value Vi at each iteration step: Vi — > Vi+i = Vi + dvi. Since Eq. 



11 ensures that both the starting vector and the source belong to ^i, v-i and dvi 
will lie also in Ai. The iteration will therefore converge to the unique solution 
within 

The solution we obtain in this manner coincides with the ghost propagator 
restricted to the subspace Ai orthogonal to the zero modes. Hence, we can 
study its behaviour in momentum space, down to (and including) the smallest 
non vanishing momentum, Pmin ~ 

2-k/L. 
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4 Results 



We calculate the ghost propagators on local sources at xq = {^,^,^,0) and 
choose the arbitrary normalization G{xq)=1 on the lattices quoted in Table |l|. 
Then we perform the Fourier transforms, for different values of lattice momenta: 

G(p)=5]G(x)e-*P-(^"-<), p^ = ^k„ k^ = 0,..,L-l. (12) 

In the practical implementation, we keep two components at /c^ = 0, and 
vary the remaining two components, covering their entire available kinemat- 
ical range. 

Figure [l| shows the resulting propagators G(p) in the case of SU(3), for the lat- 
tices 8^, 16^, 16^ X 32, and 32^. The abscissa is given in terms of the appropriate 
lattice momentum q 

q{p) = 2 ^5:sin2(^). (13) 

Some observations can be made from the ghost propagators in momentum space: 

1. On a given lattice the fluctuations of the data points appear to be small. 
In the presence of Gribov copies, however, one would expect to encounter 
large fluctuations due to gauge dependence. We therefore conclude that 
we should not worry about Gribov copies. 

2. On the different lattices all data points above pmm = 27r/L collapse nicely 
to a universal curve. Notice in particular the near-coincidence of the two 
points at p = 27r/16 which are measured in the short and long direction of 
the 16^ X 32 lattice. 

Looking at the smallest lattice 8^, we find the data point G{p = 2tt/8) to 
lie above the corresponding point from the 24^ lattice by several standard 
deviations. This tells us that finite size effects shift G upwards rather 
than downwards. 

3. Down to the smallest momentum on the 24^* lattice, p = 27r/24, the scaling 
curve increases as momentum decreases. From the non-symmetric lattice, 
16^ X 32, we can extract one data point corresponding to Pmin = 27r/32. 
This point is found to lie below its neighbours. It is unlikely that this is 
due to finite size effects as the latter would enhance G{p). 

Therefore we are tempted to interpret this point as a genuine signal, that 
G{p) has a maximum somewhere inside the interval [||^, |^]. 

It is comforting to observe a similar structure in our SU(2) data, at (3 = 2.1 (see 
Fig.^). The 16"^ and 32^ data exhibits, after due normalization in configuration 
space, small finite size effects. Again, we observe (a) agreement of G(p = 27r/16) 
as measured on the 16^ and 32^ lattices and (b) a turnover behaviour of G{p) 
at very small values of p. 



5 



We are thus led to the conclusion that there is no direct evidence for dipole 
dominance in our data. 

One may expand G{p) in powers of l/q{p)'^ 

G(p) = ^ + A + ..., (14) 

and fit the data to the lowest terms. In the region of small momenta, we were 
forced to include terms up to 0{q~^), in order to match the data. However, such 
fitting implies large cancellations of terms inducing substantial instabilities. 

In the large momentum regime , q > \/2) the data can be described by the first 



two terms in Eq.l4. The fit yields the parameters quoted in Table g. Obviously, 
the quadratic term ^ Xjq^ accounts for most of the tail in the distribution G{p), 
the effective ~ 1/?^ contribution being determined within 10% accuracy. 







a 


b 




8^ 


0.746(8) 


0.276(29) 


SU(3) 


16^ X 32(a) 


0.749(9) 


0.227(33) 




24^ 


0.744(10) 


0.256(40) 


SU(2) 


16^ 


2.073(12) 


0.267(35) 




32^ 


2.084(17) 


0.299(45) 



Table 2: Fit parameter for the q > \/2 range, according to the ansatz 14 



5 Conclusions 

We have studied, for the first time, the ghost propagator for SU(2) and SU(3) 
pure gauge theories. 

It appears that the gauge fixing algorithms used in the present study suffice 
to suppress the effects of Gribov copies. The performance of the SA algorithm 
for SU(2) is comparable to the ordinary relaxation algorithm in CPU time, but 
its has the advantage of finding the absolute minimum in a single run. The 
SOR, in combination with the Los Alamos algorithm, reduces the number of 
required iterations by nearly a factor 4, depending on w and on the lattice size. 
The extra cost for applying the "wrong" gauge transformation is negligible in 
comparison with the profit in iteration numbers. 

The ghost propagators appear to be only little affected by finite size effects as 
we can verify from comparing the data from lattices of different extensions, at 
a given /3. 

In the infrared regime - measuring down to momentum values of 0.39 GeV and 
0.86 GeV in SU(3) and SU(2), respectively - we do not confirm the expected 
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behaviour. Near the boundary of the first Brillouin zone, p vr/L, G{p) 
is dominated by 
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Figure Captions 

Figure 1: G{q) for SU(3) at /? = 6.0 on different lattice sizes. The curve 
represents the fit to the 24^-data according to the ansatz (Eq.|l4|) in the range 
q > Note that the full g-range in our simulation is given by < g < 2\/2. 

Figure 2: The same as Fig. 0, but for SU(2) at 13=2.1. The fit shown is for 
the 32^ lattice. 
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